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Abstract. Properties of nanoparticles have been studied within the framework of Ising model and the 
method of random-field interactions: the average magnetic moment and position of critical points of 
the magnetic and the concentration phase transitions depending on their size. It is shown that the Curie 
temperature is inversely proportional to the size of the particle. Critical concentration of the ordered 
state decreases with increasing of size of nanoparticles until the percolation threshold of massive 
particles. 

Introduction 

Influence of surface on the magnetic properties of nanoparticles becomes substantial with decreasing 
of their size, if the number of atoms on the surface is comparable (or exceeds) the number of atoms 
inside the particle [1]. Thus the low-temperature (97 K) measurements [2,3] revealed that the specific 
magnetic moment (//) decreases with increasing of number of atoms in nanoparticles of iron, cobalt 
and nickel up to the values of the magnetic moment of massive sample. In contrast to nanoparticles of 
3d-metals mentioned above, for rare-earth elements the magnetic moment increases with the growth 
of number of atoms [2,4]. It is shown [5], that the saturation magnetization of nanoparticles of a-Fe 
(with sizes up to 7 nm) is much smaller than in massive iron. We should note the fact, that the Curie 
temperature decreases with decreasing of the particle size [6] . 

These features of the magnetic properties of nanoparticles may be associated with a reduction 
(compared to the massive particles) in the number of nearest neighbors and, consequently, the weak- 
ening of the interaction between the atoms. 

Model 

Magnetic atoms are distributed over iV 3 sites of the cubic lattice with a probability p. Thus p is the 
concentration of magnetic atoms and N is the total number of atoms on the edge of the cubic- shaped 
nanoparticle. According to [7], the distribution function for random interaction fields H on a particle 
located at the origin can be defined as: 

W{H) = J 5 ~J2 hk ( rfc ' mfc ) II F ( mfe ) 6 ( rfc ~~ rfc °) j dr kdm k , (1) 

where 5(x — x ) — Dirac delta function, h k = h k (m fc , r fc ) — field created by atoms with magnetic 
moments m k located at coordinates r k , r k0 — the coordinates of lattice sites, F(m k ) — the distribution 
function for the magnetic moments, which in the approximation of Ising model for a ferromagnets 
can be represented as follows: 



F (mjfc) = (a k S (9 k ) + (3 k 8 (9 k - tt)) ((1 - p) 5 (m k ) + p5 (m k - m )) . (2) 



Here k — the angle between m k and OZ-axis, a k and (3 k — relative probabilities of the spin ori- 
entation along and against OZ-axis (9 k = and 9 k = n, respectively); m — the magnetic moment 
of the magnetic atom. In the approximation of nearest neighbors and the direct exchange interaction 
between magnetic atoms, the equation (1) can be represented as: 

z C z z ~ n 2 n 

W j (H) = J2p z ~ ti (1-pTY. E "iJih-M^J), (3) 

where kj is a set of nearest neighbors of the magnetic atom numbered as j, z = dim kj — its co- 
ordination number; C™ (kj) — sample of n atoms of the total number of z nearest neighbors of j th 
atom; L (Q) is a binomial set of permutations of an arbitrary set f2 with the amount of elements 
equal to 2 dimn . Introducing the notation a_ n = 1 — a n we have got = Y[lu=i-i v ekj) a ±iv an ^ 
Mi. = J2nek- i m n = m o J2nek- M^ a n — 1|- Finally, J is the constant of exchange interaction. 

Using the expression for the distribution function of the interaction fields (3), one can obtain 
equations that determine average relative magnetic moments at each lattice site: 

r ^)w 3 (H)dH = j2p z - n a-p) n i2 E ^ tanh 

B J n=0 v l v eL(C»to)) 

Expression (4) allows to investigate the dependence of the magnetic moment M = X^=i H on me 
temperature T and concentration p, as well as to determine the dependence of the phase transition 
temperature and the percolation threshold on the number N of atoms. System of equations (4) can be 
solved numerically with Newton's method. 

Results 

Fig. 1 represents the temperature dependence of the specific (per atom) magnetic moment (//) = 
M/N 3 of the particles with different number of magnetic and non-magnetic atoms. Calculations were 
performed for a simple cubic lattice. As expected, reducing the number of magnetic atoms leads to a 
reduction in the number of nearest neighbors and, consequently, to lowering the Curie temperature. 

Note that the increase in Curie temperature with increasing number of magnetic atoms is described 
fairly well by a linear function of 1/N (see Fig. 2), which is consistent with experimental results [8]. 

Dependence of the Curie temperature on the number of atoms on the edge of cubic- shaped particle 
is shown on Fig. 3. Decrease in the Curie temperature is associated both with the decrease in the 
concentration of magnetic atoms and the weakening of the interaction (decrease in the number of 
nearest magnetic neighbors). 

Reducing the concentration of magnetic atoms p and, consequently, the interaction between them 
leads to loss of ordering at a certain critical concentration p c — point of the concentration phase 
transition (see Fig. 3, on the right). One can see that the curve p c = p c (N) asymptotically tends to 
p c (N — > oo) = 0.31 — the well-known percolation threshold for a simple cubic lattice [10]. 
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Fig. 1: The dependence of the specific magnetic moment of nanoparticles of different sizes on the 
reduced temperature t = ksT '/mo J (for different concentrations p) 
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Fig. 2: Dependence of relative change in the Curie temperature on the inverse number of atoms on 
the edge 1 /N. Here Tcoo = Tc(N — > oo) is the Curie temperature of a massive particle, and ATc = 
Tcoo — Tc(N). The dots mark the experimental results [8]. We should notice that the experimental 
points correspond the lattice, that differs from simple cubic. But according to our previous work [9], 
behaviour of the curve on this figure in such the coordinates conforms to the critical exponent A and 
is independent on the type of the lattice 
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